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Abstract 

We study theoretically the formation of long- wavelength instability patterns ob- 
served at spreading of nematic droplets on liquid substrates. The role of surface- 
like elastic terms such as saddle-splay and anchoring in nematic films of submi- 
cron thickness is (re)examined by extending our previous work [Manyuhina et 
al EPL, 92, 16005 (2010)] to hybrid aligned nematics. We identify the upper 
threshold for the formation of stripes and compare our results with experimental 
observations. We find that the wavelength and the amplitude of the in-plane 
director undulations can be related to the small but finite azimuthal anchoring. 
Within a simplified model we analyse the possibility of non-planar base state 
below the Barbero-Barberi critical thickness. 

Keywords: surface-induced ordering of liquid crystals, continuum models, 
symmetry breaking transitions, anchoring phenomena, saddle-splay elasticity. 



1. Introduction 

Thin nematic film spread on liquid substrate represents a fascinating system 
to study the influence of interfaces on the organisation of liquid crystal (LC) 
molecules. The pioneering works by Lavrentovich, Sparavigna and Pergamen- 
shchik [l|, 0, 0, 3 were focused on understanding the spontaneous spreading 
of 5CB (4-pentyl-4'-cyanobiphenyl) LC on glycerol. The authors observed the 
formation of spatially-periodic stripes in films of submicron thickness resulting 
from the competition between the surface and volume energy namely between 
the boundary conditions and the elastic bulk anisotropy of nematic LC. One 
of the first models [ij, accounting for these observations, suggests the crucial 
role of the splay-bend elastic constant Ki^, which should not vary more than 
by 5 % to fit the data. Another model [3] quantifies the threshold thickness 
in agreement with experimental data, while the wavelength of stripes was not 
determined. 
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Recent systematic experiments performed by the group of Cazabat H, 0, 
H Hi suggest that the stripe phase occurs within the range of thickness ^ h ^ 
hc2, where the lower threshold hc^ ~ 20-30 nm and the upper threshold — 
0.5-0.6 iim seem to be independent of the type of substrate and LC molecules. 
The period L of the stripe domains is controlled by the thickness h of the film in 
a non-trivial way, yielding the ratio L/h oc 100 for different LC systems Q. The 
presence of the free nematic-air interface allows the thickness of the spreading 
film to be adopted by the film itself rather than being a fixed parameter. The 
striking feature of the experiments [1, H, 0, H, 0] is the existence of a forbidden 
range of thickness, visible as a discontinuity in microdroplet profile. At the room 
temperature the lower boundary is a trilayer of molecules {£.lb ~ 3.5 nm) while 
the upper boundary is (,ub ~ 20-30 nm. The structure of the film with ^ub 
is unknown, and in fact on water and glycerol striped domains coexist directly 
with trilayer and ^ub is the thickness of the thinnest striped film d, [1, Q- 
There is no theoretical approach which could account for these experimental 
observations. 

Within the continuum theory of liquid crystals [lo| nematics are charac- 
terised by the unit vector, called the director, speciiying the averaged orien- 
tation of molecules, which tend to align parallel to each other. The pres- 
ence of interfaces infiuences the surface ordering and orientation of nematic 
director, also known as anchoring, which is far from being understood com- 
pletely [ilEll- 

The conventionally assumed degenerate azimuthal anchoring 
on liquid substrates, does not necessarily mean that LC molecules can align 
along any direction simultaneously without some energy cost. Rather it means 
that because of some microscopic (e.g. adsorption layer or geometric [3] 
reasons there is a tendency of molecules to align (spontaneously or not) along 
a common preferred direction on substrate interface. To describe the planar 
degenerate anchoring at the nematic-liquid substrate interface the fourth or- 
der terms of the tensorial order parameter must be included in surface free 
energy [l5| . Although nematics are uniaxial in the bulk, the wetting layer may 
exhibit biaxiality due to the lower symmetry near the surface [l5| . According 
to [l^ the biaxiality arises naturally close to a curved surfaces with imposed 
homeotropic anchoring (free interface with air) , where the order tensor becomes 
biaxial along the principal directions of the surface. 

Our recent analysis [17j shows that without accounting for azimuthal anchor- 
ing at either of the interface, the wavelength of stripes at the lower threshold 
is infinite, in contradiction with experimental observations 0, [E, 0, The 
presence of small azimuthal anchoring at one of the interfaces penalises strong 
in-plane distortions of the nematic director (stripe domains) and favours the 
common alignment of the molecules. However, contrary to the polar anchoring, 
the effective azimuthal anchoring can be thought of as an extrinsic parameter of 
our system, induced either by the thickness gradient 1, 14, 1], biaxiality 15, 3| 
or short-range intermolecular forces. It is known [l3| for example, that part of 
nCB molecule is buried in water with a preferential orientation dictated by 
its short-range interactions with water, resulting in adsorbed monolayer, which 
then leads to the azimuthal anisotropy in nCB's orientation distribution. 
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In the present manuscript firstly we generalise the theoretical analysis ITlJlSf 
towards the distorted hybrid aligned nematic (HAN) state in presence of az- 
imuthal anchoring at nematic-liquid substrate interface. This allows us to iden- 
tify the upper threshold for stripe instability hc^ — 0.5 fj,m as well as finite 
wavelength and amplitude of the perturbation. The predictions are robust, 
compatible with experimental observations of 5CB-7CB and MBBA on glycerol 
and water and do not depend on a particular choice of elastic constants, con- 
trary to [l[ . Next we explore the implications of simplified model on the direc- 
tor's configuration in unperturbed state by i) introducing an effective azimuthal 



anchoring at both interfaces accounting for possible biaxial order [15|, [16[ or 
ii) assuming a certain scaling for the azimuthal and polar degrees of freedom of 
the nematic director. Within such an ad hoc approach we find twisted state and 
modification of the Barbero-Barberi condition [3]. Clearly these approaches 
are oversimplified to understand the physical picture as a whole, where the in- 
terplay between hydrodynamic and elastic instabilitites [2^, [U , may lead 
to a nontrivial film profile extending over the micrometers and thickness dis- 
continuity, which occurs within nanometers scale @, Nevertheless, it 
is a good starting point to investigate the connections between molecular and 
mesoscopic scale phenomena quantitatively. 



2. Formulation of the problem 

Nematic liquid crystals are described by the unit vector n, called the director, 
characterising the averaged orientation of molecules. The free energy associated 
with distortions of this director in space is given by 

•^e;-^ / dF{i^|Vn|2-X24V-(n(V-n)+nx Vxn)}, (1) 

where we have assumed the one-constant approximation of the Oseen-Zocher- 
Frank free energy flo|. quadratic in the director derivatives. The elastic constant 
K stands for the equal splay, twist and bend elastic moduli [23], while the saddle- 
splay elastic constant if 24 should satisfy the Ericksen's inequalities 24 1 



\K24\ ^K, K^O, (2) 

which guarantee the stability of the uniform nematic state n ~ const. The 
uniformly aligned director configuration can become unstable in presence of the 
competing boundary conditions and the saddle-splay elastic term, favouring the 
distorted configuration of the nematic director. Therefore the modulated stripe 
state can be achieved, when the saddle-splay contribution to the free energy 
becomes substantial, which is plausibly realised in thin nematic films with weak 
anchoring boundary conditions. 

For convenience we parametrise the director as 

n = sin 9 cos + sin 9 sin (pey + cos 9 e^ , (3) 
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where 6 and are the polar and the azimuthal angles, depending on the Carte- 
sian coordinates x,y,z. The first term of the free energy ([T|) integrated over the 
volume V characterises the bulk contribution and within the parametrisation 
([3]) can be rewritten as 

^6 = y^d^(|V0p+sin20|V(pp), (4) 

while the last divergence term can be transformed into the surface integral as 

Ts — / rfS* [cos</j(92:6' + cos0sin6'9j,(/3) + sin(p((9y6' — cos0sin6'9j;(^)] 

(5) 

where h is the thickness of the nematic film and the surface normal in fiat central 
part of the spread film is assumed to be parallel to the z-axis. 

We assume that the anchoring energy admits the general Rapini-Papoular 



form [25 1 



^a = \Y^ [ dS{ Wg^ sin^ {9, -e,) + W^^ sin2 9, sin^ - ) } , (6) 

1=1,2 •^•S' 



where the integration is done over the nematic liquid substrate interface i — 1 
{z = 0) and the nematic-air interface i — 2 {z — h). Here Wg^ is the po- 
lar anchoring strength, accounting for energy cost due to the deviation of the 
polar angle 6i {9i = 9{0) and 92 = 0(h)) from the easy axis 9i = 7r/2 (pla- 
nar anchoring) and ^2=0 (homeotropic anchoring), respectively. The second 
term in J'a accounts for the deviation of the director's projection on the xy- 
plane from a common alignment, which happens during the formation of stripes. 
From the experimental works [j, '^6 know that the planar anchoring Wg-^ is 
stronger than the homeotropic anchoring Wg^, and quantitatively the extrap- 
olation lengths for 5CB on glycerol system are Lg^ — K/Wg^ — 0.7 fim and 



Lg-^ — K/Wg^ — 0.35 /im. According to the value of the azimuthal an- 
choring Wif,-^ 2 should be one or two orders of magnitude smaller than the polar 
anchoring, but in general it is not fixed a priori. 

In the following sections we assume the vanishing azimuthal anchoring at the 
nematic-air interface W^p^ = 0, while W^p^ ^ 0, and explore the predictions of 
the formulated model from the bifurcation analysis point of view. In section [??T] 
we extend theoretical analysis towards the alternative case, when W^.^ ^ and 
indicate the consequences of the proposed model for the director's configuration, 
without referring to any particular experiment. 



3. Linear stability analysis 

In this section we assume the vanishing azimuthal anchoring at the nematic- 
air interface W^p^ = 0, yielding base states with ip = (pi = const, confined to 
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planar hc^ HAN h 




Figure 1: Schematic representation of the structural phase transitions between the planar 
state, striped state (/ic^ < h < h^^) with periodicity L in the j/-direction, and the hybrid 
aligned nematic, where the film thickness h is thought of as an order parameter. The critical 
thickness between the planar and HAN state is hc^ = Lg^ — Lg^ [19| . 



the 2D plane. Hence, the angle 6{z) depends only on the coordinate along the 
thickness of the film z — z/h and satisfies the Euler-Lagrange equation 

fe0 = O, yielding e{z) = 0^ + {62 - 9^)2. (7) 

The equilibrium angles 9i, between the director n and the z-axis at liquid sub- 
strate z = and at air ^ = 1 interfaces should satisfy the following boundary 
conditions for a given thickness h 

2(6li - 6I2) - sin 26*1 = 0, (8a) 
Let 

2(01 - 6*2) - sin 26*2 = 0. (8b) 

We will distinguish two base states 

• the homogeneous planar state with 9 = 9i = 6*2 = 7r/2 {h < hcg), 

• the HAN state with 0i 02 {h > hc„) satisfying ^ with 6 

The critical thickness hcg = Lg^ — Lg-^ (Barbero-Barberi [l^) characterises 
the transition between the planar and the HAN states. To find the critical 
thresholds for the formation of stripes we consider the variation of the total 
free energy T — Fb + Ts + J-a given by the equations (U])-® with respect to 
the small perturbation of the director n ([3]) . The perturbation angles -0 (polar) 
and (azimuthal) are assumed to be small 0{e) with periodic modulation in 
2/-direction 

'<P{y. z) = ef{z) sm{xy), (j){y, z) = eg{z) cos{xy), (9) 
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where y — y/h is the dimensionless coordinate, and x ~ ^irh/L is the wave- 
number, with L being the period of stripes. Close to the instabihty threshold 
we can introduce a small parameter e, defined as h = /ic(l +e^), where he is the 
critical thickness, characterising the instability of planar (HAN) states towards 
the stripe state. The difference in the total free energy of the stripe state and 
the base state can be written as the Taylor series 5T = e^ST^^'' + e^ST^'^^ + . . ., 
with 



dy 



dS { sin2 0{S)[{d,cj,)^ + {dycj^f] + {d,4'f + {dy4'? } + 
h cos 292 



L92 



-TP{iy + 21} sin^ dy(j> 



(10) 



where i9 = K2i/K is another dimensionless parameter related to p = {K22 + 
K24:) /K in [l^ as I? — p—l. The value of ^^24 as well as the extrapolation length 
for azimuthal anchoring Lip^ are not precisely known, therefore we present below 
our results for the critical threshold as function of i? and L„ 



The Euler-Lagrange equations associated with (ITOl) are 

d,Bf - xV = 0, 
aj(sin2 ed^g) -x^sinHg^O, 

with the general solution given by: 

f{S)=Cie^' + C2e-^', 

9V'> = ■ at~\ ' 



(11a) 
(lib) 



(12a) 
(12b) 



where 9{z) = 6'i + {62 — 6*1)2 and A = \Jx^ — [62 — OiY- In the planar case 
6 ~ t: /2 equation (|12pb simplifies to the sum of two exponents. The integration 
constants d can be found from the boundary conditions associated with ([T0| 



(2) 



(2) 



(2) 



doj. 



(2) 



9V 



= 0, 



0, 



duj, 



(2) 



dy 

d dcji"^^ 
dy d(t)^y 



(2) 



= 0, 
(13a) 

0, (13b) 



6 



which can be cast into the matrix form as X]i=i -MijCi = 0, where 
/i cos 2^1 hcos2ei i?;^sin6'i 



M - 



q h — L^^ A4--£/^i {02-~0i) cot Qi h-\-L^^ A+L^^ {O2—O1) cot ^1 

sin 9\L^p-^ sin 6\L:^-^ 

_a^-X e^(A-(i92-6'i)cot 6>2) e"^ (A+C^s -6*1) cot ^2) 

sin 02 sin ^2 

(14) 




Figure 2: The critical thickness, characterising the lower (/ici) and the upper (hc^) instability 
thresholds as function of saddle— splay elastic constant i? = K2A/K. The stripe phase exists 
between the curves of the same colour and it is restrained when L^-^ decreases or azimuthal 
anchoring at nematic-liquid substrate interface W,p-^ increases. The values for the polar 
anchoring are known Lg^ = 0.35 piva and Lg^ = 0.7 pLxa f^. |26|| . 

The nontrivial solution (C; 7^ 0) exists if and only if the determinant of the 
matrix is zero, det A4 = 0. This condition gives the relationship between 
the thickness h of the nematic film, the dimensionless wavenumber the angles 
01^2 (HAN case) and other parameters of the system. Since the critical thickness 
of the nematic film he is unknown, the problem should be solved consistently 
together with the boundary conditions for the equilibrium angles 9i ([8]), which 
in turn also depend on h. The minimum of the curves det = in the h-x 
plane allows to identify the critical thickness he and the critical wavenumber Xc 
at the bifurcation point, for a given value of i?, and the extrapolation length of 
azimuthal ancho ring L^^, while the polar anchoring is known Lgi = 0.35 /zm, 
Lgj = 0.7 /zm In Fig. [2] we plot h^ characterising the planar-stripe (/ici) 

and the HAN-stripe (/1C2) instability thresholds. The curves are symmetric 
around "d — Q. In presence of strong azimuthal anchoring (small values of 
L^J the formation of a stripe phase is suppressed, as was predicted in (27j . 
Experimentally compatible values for critical thickness /ici — 40 nm and — 



7 




Figure 3: The critical wavenumber Xci and Xc2 > characterising the lower and the upper 
instability thresholds, respectively. In the limit of the vanishing azimuthal anchoring — >■ 
oo the wavelength of stripes vanishes. 



0.5 fim occur at ~ 1 and for of order of tens of micron, corresponding 
to extremely weak azimuthal anchoring. 

In Fig.|3]we plot the critical wavenumber Xci at ha, with one inverse branch 
for illustrative purposes. As mentioned above, the experimentally relevant crit- 
ical thickness ha occurs at ^^24 — K and L^-^ ~ 70 /zm, yielding the critical 
wavenumbers Xci — 0.2 and Xc2 — O-l- For planar-stripe instability x vanishes 
only in the limit of — ?> oo {Wip-^ — 0), while for HAN-stripe instability x — ^ 
is always the solution of det^Vl = 0. Comparing the results on h^ and Xci with 
our previous analysis [l7j . where we have taken into account the twist elastic 
constant if 22 7^ ^11 and assumed a non-zero azimuthal anchoring at nematic- 
air interface, we conclude that the predictions of both critical parameters agree 
quantitatively. This plausibly means that the ratio K22/K11 does not play a 
dominant role in our system and that the one constant approximation is reason- 
able to analyse the stripe instability. It is hard to directly compare theoretical 
predictions of the upper threshold hc2 and Xca with experimental data, because 
of the presence of defects in nematic films thicker than h ~ 0.5 /Ltm 

UMM- The 

HAN state, therefore, is an idealised director's configuration which is probably 
not accessible in real systems, because of the presence of defects. For complete- 
ness in Fig.|4]we plot the equilibrium angles 61^2 found at hc2 and satisfying the 
boundary conditions ([5]). 

4. Weakly nonlinear analysis 

The linear stability analysis, considered above allows to identify the critical 
thickness and the critical wavenumber at the instability threshold. However 
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Figure 4: The equilibrium angles Si, 2, satisfying the boundary conditions l(8} (for = 0) at the 
critical thickness fee, (see Fig.[2]l. The values for the polar anchoring are known Lg-^ = 0.35 
and Lg^ = 0.7 ^tm 0,1111 • 



the linearised boundary conditions (|14|) give only the ratio of the integration 
constants d. Thus, the perturbation functions and ip are known up to some 
multiplying factor A, called the amplitude of perturbation. To eliminate this 
uncertainty and to find A we expand the total free energy further in powers of e 
up to O(e^), fourth order in (j) and ip. Moreover, we assume to be in the vicinity 
of the threshold he, Xc, and the thickness of the film given by ft. = hc{l + e^) is 
the control parameter of our system, with e ^ 1 measuring the distance from 
the threshold. Then the Taylor expansion of the free energy is 

where ^J^'^^^ is given by (flU)) for planar (HAN) case, dJ-a'^'' is the second order 
contribution of the anchoring energy and ST^^"^ is the fourth order contribution 
of the total energy, expanded in (j) and ip. We make use of Mathematica to 
integrate directly (fT^ . since the functions (j) and as well as the critical 

point were found before in section [3l The first term vanishes at the critical 
point SJ^^^^ \h=hc, = 0, as expected. Then the amplitude can be identified by ex- 



tremising the next order term, namely = —bT^a^j (2(5 J^*^^^) . As pointed 

I , ^ ^ _ _ X=Xc 

out in 28] this procedure may be equivalent to the derivation of the amplitude 
equations. The energy method, however, does not allow to find a nonlinear 
change in the wavelength. 

In Fig. [S] we plot the amplitude A for the planar-stripe instability as func- 
tion of the azimuthal anchoring L^p^. The increase of A with L^^ suggests 
the suppression of the stripe phase in presence of strong azimuthal anchoring, 
which also agrees with the linear stability analysis performed in section [3l For 
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|A| 




Figure 5: The amplitude for the planar— stripe instability as function of the azimuthal anchor- 
ing -Lipi- The choice of parameters is Lg-^ = 0.35 fira, Lg^ = 0.7 fim, = 0.95, = and 
the corresponding thresholds are found in section [3] 



all considered values of L^^ the amplitude is 0(1), justifying the perturbation 
approach developed above. 

The analysis developed above predicts the threshold for stripe instability, 
compatible with experimental observations, the question about the base state, 
however, remains a sensitive issue of the model as mentioned above in the in- 
troduction. The assumption about 'idealised' base states, such as planar and 
HAN, is not strictly supported by recent experiments, since the structure of the 
film with 20-30 nm thickness is not known while at the upper threshold « 0.5- 
0.6 /im stripes are replaced by other complex structures with defects 0, H, Q. 
The alternative 2D model for the spreading of nematic drops accounting for 
defects and the evolution of profile was proposed in 2l|, |22| , where an ad hoc 
anchoring condition which relaxes anchoring strength when thickness of the film 
goes to zero was introduced. Indeed, the Rapini-Papoular anchoring condition 
is plausibly not suitable for an ultrathin nematic films with a free interface. 



5. On the base state of nematic films 

In this section, without referring to any particular experiment, we explore 
the possibility of different director configurations along the thickness of the film 

i) by imposing an effective azimuthal anchoring at both interfaces ([6|) (Wip^ ^ 0), 

ii) by assuming a certain scaling of the in-plane and the out-of-plane degrees of 
freedom, leading to their coupling. Then within the developed ansatz we look 
for the critical thickness between the planar state and distorted state, yielding 
the modified Barbero-Barberi condition hc^ = Lg^ — Lg^ jl9l] . 

5.1. In presence of azimuthal anchoring W^^ 

The presence of effective azimuthal anchoring might account for biaxiality of 
nematic order close to the surface 15, 3l, thickness gradient [l,[l3,B or short- 



range interactions [13|, |12| . We make use of the continuous approach formulated 
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above in the section[31 For flat films the angles 9{z) and (p{S) ^ depend only on 
the coordinate z = z/h along the thickness of the film and the Euler-Lagrange 
equations associated with the bulk free energy ^ are 

fe6i -sine cos 6l(9j(^)2 =0, (16a) 
di{sin^ ediip) = 0. (16b) 

The first integral is 

9^^ = -%, (9,0)2 = s|__4-, (17) 
sm sm 

and the second one yields the solutions and ip, written in the following form 



cose{S) = -^^-^ L sin(B3 - B2S), tan{ip{z) + B4) = — i tan(S2Z - B3). 

B2 B2 

(18) 

Let us consider several limiting cases: 

a) 61 — 02 — tt/2 and (pi — (p2, (Si = 0, B2 = 0), planar state with uni- 
formly aligned director n, 

b) e{z) = 01 + {82 - 0i)z and (fi = Lp2, {Bi = 0), HAN state with varying 
along the thickness of the film 

c) 01 = 02 = 7r/2 and (p{z) = ipi + {ip2 — >fi)z, (i3| — Bf), twisted state 
with If varying along the thickness, 

d) 01 ^ 02 and fi ^ <^2, satisfying ([TSl) . (fT9|) . twisted— bent. 

In general, the integration constants Bi should be determined from the 
boundary conditions associated with Qj, ^ expressed as 

2526*1 _ + sin 20i(— — sin^ ipi] ^0, (19a) 

h h 



2di9\^^^^+sm202{— + —8m'{ip2->p)] =0, (19b) 
sin2 0i(^2c>i¥'|j^o - sin 2^1^ = 0, (19c) 
sin2 6i2(^295V7|j^^ + -^sin2 6l2sin2((^2 =0, (19d) 

where Lg. and L^. are the extrapolation lengths for polar and azimuthal an- 
choring, respectively, Lg^ > Lg-^ (the planar anchoring is stronger than the 
homeotropic one), as before we assume 0i = 7r/2, 02 = and the system of 
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coordinates is chosen so that </?i = 0, and then the index is dropped by for 

(pi = (p. 

The conventional uniform planar and HAN ground states can be found for a 
vanishingly small azimuthal anchoring L^^^ ^ Lq^ and <^ = 0, 7r/2. Interestingly, 
if e (0, 7r/2) is an arbitrary angle, the boundary conditions are satisfied if 
and only if ipx =^ ^2 or d\= Qi = 0, irrespective of the smallness of the anchoring 
strength W^^. The planar state becomes linearly unstable with respect to the 
twisted state if L^^ — L^^ < h < Lq^ — Lg^. The states a) and c) become 
unstable with respect to the small perturbations of the polar angle 6, resulting 
in states b) and d), respectively, if the thickness of the nematic film satisfies the 
following inequality 

/i(L^j - sin^ + Lg.^ sin^((p2 - (p)) > L^^L^^ x 

(Lea - Loi) - Le^Lg^{L^p^ sit? (pi + L^-^ sin^((^2 - <^))- (20) 

In the limit of the vanishing azimuthal anchoring (|20p yields the known result for 
the critical thickness hcg = Lg^ — Lg^ [l^ . In presence of azimuthal anchoring 
the threshold /ico cannot be defined explicitly, because of the dependence of i^i 
and ip2 on h through the boundary conditions (|T9l) . In Fig. |6]we plot the com- 
puted threshold ft-coj which has a smaller value than the Barbero-Barberi result, 
meaning that in presence of azimuthal anchoring the homogeneous planar (or 
twisted) states are easily replaced by the HAN (or twisted-bent). Experimen- 
tally, he structure of the film with thickness 20-30 nm is unknown Q and 
one cannot exclude the tendency of the director to form a twist configuration. 
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Figure 7: The relative stability of the planar state w.r.t. perturbed state lj22p {Lg-^ = 0.35 fim, 
Lg^ = 0.7 iim) a) e^-term in II23I I is always negative above the critical value of the in-plane 
95-undulation rjc = arcsin[(Lg2 — Lg^ )/{Lg^ +Lg^ )]; b) under the bell-shaped curve < I I23I I 
the planar state is linearly unstable. We choose e-^ = 0.5, for smaller e the peak of the curve 
shifts towards smaller r] and the area under the curve shrinks. 



5.2. Different scaling for 9 and Lp 

Thin films with planar configurations 9 = tt/2 and = 0, can become 
unstable with respect to a small perturbations 9 and (p, admitting the following 
form 

9(i) = ^+e9 + ^{z) = ip + e^ipi + ... (21) 

with £ being also small. Hence we have conjectured that the reorientation of n 
in the plane of the film happens before the change of the polar angle 9. Indeeed, 
according to 11] as well as our analysis, the azimuthal anchoring should be 
one or two orders of magnitude smaller, resulting in energetically less expensive 
reorientation of the molecules in the plane of nematic-liquid substrate interface 
rather than change of the polar angle 9. Usually, (p is considered to be negligible, 
resulting in the planar ground state and thus impossible coupling of two degrees 
of freedom. In our case, the solution of ([T6|) up to 0(e) takes the form 



(p{z) ^ri- z, e{z) = A cos{r] ■ z) + B sm{r] ■ z), (22) 

with coupling of and 6 resulting in negative contribution to the linearised 
total free energy 



2h 



v'+e'{ I dzm9f~ri'9^} + ^9\Q)-^e\l) 







(23) 



Although the first term is always positive the second term becomes negative if 
the thickness of the film is ft, < /i_ oi h > h^, where 



2 sm rj \ * / 
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yielding the Barbero-Barberi critical thickness hc„ — Lg^ — Lg^ [l^ in the 
hmit of 77 — >■ 0. The curve h±{rj), shown in Fig. [7^, separates the region of 
the unperturbed planar state from the region where the instabihty towards the 
state with non-zero (p and 9 may plausibly occur. To quantify the destabihsing 
effect of the second term compared to the stabihsing effect of the first term we 
find the amplitude of the perturbation as described in the previous section |31 
In Fig. we plot the curve when = so two competing effects cancel each 
other, giving the bell-shaped region for the planar state instability. 



6. Discussion and concluding remarks 

One of the main challenges in modelling thin nematic films is to describe the 
interplay of different length scales simultaneously Q . The profile of nematic 
microdrop extends over micrometers while the thickness of the film experiences 
sudden jump from 3 nm to 30 nm. Similarly, the wavelength of observed stripes 
varies from microns till hundreds of microns while the thickness of the flat 
central part of nematic film changes in the range 20 nm < h < 0.5 /im. For such 
small thickness the interfaces play a crucial role, and short-range interactions 
should be included in theoretical models. This would require the development 



of microscopic approaches, probably similar to 1^, which would result in better 
understanding of wetting and anchoring properties of nematic liquid crystals. 
On the other side of the length scale, it is useful to account for the film profile by 
studying hydrodynamic (in)stabilities as described in 2^, 22 1. Phenomenology 



is an alternative approach which can be useful to capture both long-wavelength 
instability patterns as well as the structure of the ultrathin nematic film and to 
interprete experimental observations. 

Here, within a phenomenological framework, we have reconsidered the stripe 
instability in thin nematic films, subjected to antagonistic boundary conditions. 
The thickness of the film as well as azimuthal anchoring were treated as ex- 
trinsic parameters of our system related to a particular experiment, while the 
polar anchoring and the saddle-splay elasticity being intrinsic ones. Solving the 
variational problem within the one-constant approximation we identified the 
critical thickness and the critical wavenumber, which determine the threshold 
of stripe instability. We (re)examined the effects of saddle-splay constant K24 
as well as azimuthal anchoring on the critical thickness, wavelength and the 
amplitude of stripes at the bifurcation point. Then, we analysed the possibility 
of non-planar base state below the Barbero-Barberi critical thickness, assuming 
a certain anzats for the biaxial order of nematic director in the vicinity of the 
interface or a certain scaling for the in-plane and out-of-plane degrees of free- 
dom. Considering different base states and their perturbations is relevant for 
predicting instability patterns such as stripes, squares and chevrons as well as 
for establishing further connections with experiments. 
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